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Abstract 

y—{ Charged B decays to three charged kaons are analysed in the framework of the QCD factorization approach. The strong final state 
7— I K + K~ interactions are described using the kaon scalar and vector form factors. The scalar non-strange and strange form factors at 
low K + K~ effective masses are constrained by chiral perturbation theory and satisfy the two-body unitarity conditions. The latter 
stem from the properties of the meson-meson amplitudes which describe all possible S -wave transitions between three coupled 
>~>channels consisting of two kaons, two pions and four pions. The vector form factors are fitted to the data on the electromagnetic 
^ kaon interactions. The model results are compared with the Belle and BaBar data. Away from 0(1020) resonance, in the S -wave 
^> dominated K + K~ mass spectra, a possibility for a large CP asymmetry is identified. 
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1. Introduction 

Recently, charmless three-body decays of B mesons have 
i been intensively studied both experimentally and theoretically. 
On the experimental side, Dalitz plot analyses of the charged 
B decays were performed by Belle U] and BaBar Q collab- 
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orations. Likewise, several theoretical studies involving the 
B ± K+K-K* decays have been published |H, H and fl. 

Since charged kaons interact strongly, their long distance in- 
teractions in the final states have to be well understood if one 
aims at extracting weak decay amplitudes from the B to KKK 
decays. In this Letter we go beyond an isobar model parameter- 
ization of the B decay amplitudes and introduce additional the- 
oretical constraints on the S -wave two-body K + K~ interaction 
amplitudes, which follow, in particular, from unitarity. In order 
to satisfy unitarity in two-body interactions we construct scalar 
strange and non-strange form factors which enter into the S- 
wave parts of the decay amplitudes. These amplitudes are cal- 
culated in the framework of the QCD factorization approach. In 
the construction of form factors we use experimental informa- 
tion on the K + K~ interactions coming from experiments other 
than B decays, for example from K + K production processes 
in hadronic collisions or from <? + e~ reactions. We apply also 
some low-energy constraints coming from the chiral perturba- 
tion theory. Preliminary results of our analysis concerning the 
B ± — > K+K K* reactions can be found in Ref. |0]. 

In Section 2 we formulate the theoretical model of the B + and 
B decay amplitudes. Presentation of results and their compar- 
ison with the experimental data are given in Sec. 3. Our con- 
clusions are presented in Sec. 4. 
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2. B ± -> K+K-K* decay amplitudes 

Inspection of the Dalitz plots of the Belle IH and BaBar Q\ 
experiments reveals an accumulation of events for the K + K~ 
effective masses below 1.8 GeV. Indeed, several mesonic reso- 
nances which can decay into the K + K~ pairs exist in this range 
J3l. Among them there are scalar and vector resonances which 
are formed via the S - and f-wave final state interactions. In the 
first approximation one can neglect their interaction with the 
third kaon. This justifies using the QCD quasi-two-body fac- 
torization approach for the limited range of the effective K + K~ 
masses (see, for example Ref. [8]). The B — > K + K~K~ ampli- 
tude is then expressed in terms of the following matrix element 
of the weak effective Hamiltonian H: 

{K-( Pl )K + (p2)K-(p 3 )\H\B-) = A s + Ap, (1) 

where the S -wave part is 

G F < fl 



(m z K )yr n 2 (s 23 ) 



25, 



(2) 



+ ^^(M 2 B - m 2 K )F$ K ( S23 )vr((s 23 )}, 
m b -m s > 



the F-wave part is 



A P = ^L{^A B p (ml)yFf K -(s 23 ) - F»\ S23 )\ Wu F? K \s 23 ) 
V2 Jp 



+ w d Ff K ~(s 23 ) + w s Ff + *~0s 23 )]}4~?i -~p 2 



(3) 



and the interacting kaons are taken to be kaons 2 and 3. Fur- 
thermore, s 2 3 is the square of the K + (p2)K~(pi) effective mass 
m23 = mic+K-, while "pi and~/?2 are the kaon 1 and kaon 2 mo- 
menta in the center of mass system of the kaons 2 and 3. The 
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scalar product of the kaon momenta can be written in terms of 
the helicity angle 



~Pl -~P2 = -\~Pl\\~P2\C0S® H . 



(4) 



In these equations Gf is the Fermi coupling constant, = 
0.1555 GeV and f p = 0.220 GeV are the kaon and the p me- 
son decay constants, Mp, Mk, nib - 4.9 GeV, m s = 0.1 GeV, 
m u = 0.004 GeV and ma = 0.004 GeV are the masses of the 
B meson, kaon, fe-quark, strange quark, down- and up-quarks, 
respectively. 

The functions F" and Tl, present in the S -wave amplitude 
in Eq. (0, are the kaon non-strange and strange scalar form 
factors. The vector form factors F^ (for q — u,d and s), 
introduced in Eq. ©, are defined through matrix elements 



In terms of the quantities introduced above one defines: 

y = A u [a Xy + a\ y + a" Wy - (a u 6y + < y )rf ]+ 



where 



A c [a c 4y + a\ ()y - (a c 6y + cfyr*], 



x (m h + m u )(m u + m s ) ' 



(6) 



(7) 



w u = AJa 2w + a iw + a 5w + a 7w + a 9w ) + A c (a iw + a 5w + a 7w + a 9w ), 

(8) 



Wd 



= A M [a 3w +a5 Vl ,-i(a7 W -i-a9 W )]-i-A c .[a3 W -i-a5 W --(a7 W -i-a9 W )], 

(9) 



< K + (p 2 )K-(p 3 )\q 7tl q\0 >= (p 2 - P3 ) ll Ff K '(S23), (5) 

where |0 > is the vacuum state. The K + K pair in the S -wave 
is then denoted by R s = (K + K-) s . Similarly /? P = (K + K-) P 
stands for the P-state. Furthermore F^^ K K >s in Eq. d2j is the 
form factor of the transition from the B meson to the K + K~ 
pair in the S -state, % is the constant related to the decay of 
the (K + K~)s state into two kaons, and Bq = m^/(m u + mj), 
where m n is the pion mass. We take F^T^ K K ^ s (m 2 K ) = 0.13 
|@] and we fit x to the data. Functions Fq K (s2i) and Ff K (s23) 
are the B -> K scalar and vector transition form factors and 
Aq P (;«^) = 0.37 H is the B — > p transition form factor. In our 
approximation, the ratio A^ p /f p represents a general factor re- 
lated to the transition from B to any (K + K~)p state and then its 
decay into the final K + K~ pair. For the case of the p meson this 
coupling to the pair of kaons is effectively realized only above 
the K + K threshold. 

The weak decay amplitudes depend on QCD factorization 
coefficients a p - and on the products A„ = V u bV* s , A c = V c bV* s , 
where Vij are the CKM quark-mixing matrix elements. In or- 
der to describe B decay into mesons M\ and M2 we follow 
Ref. and calculate the coefficients a p .(M\ M2) at the next- 
to-leading order in the strong coupling constant at the renoma- 
lization scale equal to nib/2. Here the M\ meson has a com- 
mon spectator quark with the decaying B meson. In the case 
of the B — > K + K K decays, M\ or M2 can be either kaon 
K~, or systems Rs, Rp. We take into account one-loop vertex 
and penguin corrections to a p {M\M2) but neglect those due to 
hard scattering or the annihilation since they are expected to 
be generally suppressed. In the QCD factorization approach 
they receive logarithmically divergent contributions due to soft 
gluon interaction which are "unavoidably model dependent" 
(see Ref. 18D). We treat such soft interactions by introducing the 
form factors constrained by data on meson-meson interactions, 
taken from analyses of reactions other than the B decays. Un- 
der these conditions we have a P j(RsM2) = a p .{RpM2), with their 
common value denoted below by a p (Rs,pM2) = a p y . We also 
use the abbreviations: a p = a p (K~Rp) and a p = a p (K~Rs). 
The values of coefficients a p (M\M2) are given in TableQ] 



A„[a 3w + a\ w + a 5w - ^(a 7w + a % , + a" 0R ,)]+ 
A c [a 3w + a c iw + a 5w - ^(a lw + a 9w + a c Ww )j, 



(10) 



and 



V = A u (-a u 6v + l -al) + A c (-al + (11) 

One can notice that in the expressions for the decay ampli- 
tudes there are no transitions to the K + K~ states of spin 2 or 
higher. This results from the application of the factorization 
approach in which matrix elements to spin states higher than 
one vanish. The contribution of /2(1270) with its rather small 
branching fraction to KK (4.6 %) is thus not included in this 
study. 

Since two identical charged kaons appear in the final state of 
the B — » K + K~K~ decay, the amplitude of Eq. ([TJ has to be 
symmetrized 

A- ~\{K-( Pl )K + (p2)K-(p 3 )\H\B-)+ 

V2 L (12) 



{K-(p 3 )K + ( P 2)K-( Pl )\H\B-)]. 



The symmetrized amplitude for the B + — > K + K K + reaction 
reads 



Km = A sym {A u -> A*„ A c -> A^, B 



B + ). 



(13) 



The final state kaon-kaon S -wave interactions are dynami- 
cally coupled with systems consisting of two and four pions. 
Thus a system of three coupled channels: nn, KK and An (effec- 
tive (27r)(27r) or crcr, pp etc.), labelled by j = 1,2, 3, is consid- 
ered in the construction of scalar form factors F!J and Tl. Here 
we use an approach initiated in yj and recently developed in 
ifioll for the B ± — > Ti^iCn* decays. A set of the 3x3 transi- 
tion amplitudes T, describing all possible transitions between 
the three channels, is taken from a unitary model of Ref. Ill 111 
(solution A). We introduce two kinds of production functions 
R": s , labeled by n (non-strange) or by s (strange): 



R"j' s (E) 



+ t"' s E + oj": s E 2 
j j 

1+cE 4 



1,2,3, 



(14) 
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where a"' s , t"' s , ofj s and c are constant parameters, while E rep- 
resents the total energy and is related to the center of mass mo- 
menta kj = ^E 2 — m 2 , with m\ = m n , n%2 — trig, mi, = 700 

MeV, and s = E 2 = m 2 K+K _. The three scalar form factors, writ- 
ten in the compact row matrix form Y n ' s *, are given by 



r- s * = R n ' s + TGR" 



(15) 



where R n,s are rows of the production functions and G is the 
matrix of the Green's functions multiplied by the convergence 
factors Fj(p) = (k 2 + k 2 ) / (p 2 + k 2 ). These factors, which reduce 
to unity on shell (p = kj), make finite the relevant integrals over 
the intermediate momenta p. The parameter k will be fitted to 
the data of the BaBar [2] and Belle [1] Collaborations. 

For both the non-strange and strange form factors we also 
constrain their low energy behaviour using the chiral perturba- 
tion model of Refs. 1 12, 13ll . At low s values one writes the 
following expansion: 



T n '\s) s d": s + 
j K ' j j j 



s, y = l,2,3, 



(16) 



with real coefficients d"j s and fj' s . Explicit formulae for the set 
of non-strange form factors, in particular for the presented in 
Eq. (0, are given in Eqs. (24-35) of Ref. fioll . For the strange 
form factors we have 



dl ~ 2 



h - 2 



16m: 



f 

8^ 

f 2 



72tt 2 / 2 



„2 \ 



1+log-f 



32n 2 f 2 



.2 \ 



1+log^ 



432n 2 m 2 f 2 



(17) 



(18) 



and 



4 = 1 + 



8(2^ - L' 4 ) 



32L, 



f 2 



8 2 



4%n 2 f- 



(ml+4m 2 K )- 



log^ + 



i6l: 



5 2 



36tt 2 / 2 



,2\ 



1+log^ 
fi 2 



(19) 



4l: 



h - f 2 f 2 



2\6n 2 f 2 m 2 32n 2 f 2 



,2\ 



1 + log -i 

V 2 



64n 2 f 2 



,2 \ 



(20) 



1+log^- 

p 2 



In these equations m n is the 77 meson mass, p. is the scale 
of the dimensional regularization and / = f n / y2. Using 
/ = 92.4 MeV and the chiral perturbation theory constants L r k , 
k - 4,5, 6, 8, given in Table X of Ref. 11411 . we obtain the non- 
strange-sector parameters: d" = 1.1957, /" = 3.1329 GeV 2 , 
d\ = 0.7193 and j% = 1.6719 GeV* 2 and their strange-sector 
counterparts: d[ = -0.0016, // = 0.2393 GeV 2 , d\ = 1.0410 
and / 2 ! = 0.6235 GeV 2 . For the form factors related to the 
third channel at low energies we make the simplest assumptions 
d'l = d-; = f» = // = 0, as in Ref. HI. 

The coefficients a"j'\ t".' s and u>"j s are constrained by the val- 
ues of the form factors at low energies. They are calculated 



using the low energy expansion of Eq. (15[ and are listed in 
Table [2] The parameter c, which controls the high energy be- 
haviour of R, is fixed while fitting the data. 

Our scalar form factors satisfy the following unitarity condi- 
tions: 

Imr* = T'Dr*, (21) 

where D is the diagonal matrix of the kinematical coefficients 
which are proportional to the channel momenta kj in the center 
of mass frame: 



Di 



kys/s 

8n 



6ij8(?Js-2mj), i,j= 1,2,3. (22) 



Presence of the resonances in the K + K effective mass distri- 
butions (see Refs. [ 2L LjJ] ) is a direct manifestation of the K + K~ 
final state interactions. The most prominent resonance in the P- 
wave is 0(1020). In 2005 Bruch, Khodjamirian and Kiihn [15] 
described the electromagnetic form factors for charged and neu- 
tral kaons in terms of additive contributions from eight vector 
mesons: p = p(770), p = p(1450), p" = p(1700), u = w(782), 
o) = w(1420), lo" = w(1650), <p = 0(1020) and <p' = 0(1680). 
Using quark model assumptions and isospin symmetry as in 
Ref. [15] one can deduce the following expressions for the three 
f-wave form factors F% K defined in Eq. (0: 



7 K*K- 



1 



7 K + K- 



= -(c p BW p + c pl BW p ,+c p „BW p ,> 
+ Cld BW w + cjBWj + cj'BWj'X 



= -(-c p BW p - c-BW > - c "BW » 

2 x y y p p p p (24) 

+ c u B? w + c,j BW,j + c,j- BW,,/- ), 



(25) 



In the above equations BWj, i=l,...,8, are the energy-dependent 
Breit-Wigner functions, defined for each resonance of mass m, 
and width T, as 



BWi(s) 



1 2 - s - i -\fs F;(.s) 



(26) 



and c, are the constants given in Table 2 of Ref. 11511 for the 
constrained fit. 

The B to K transition form factors have been parametrized 
according to Ref. 01611 : 



F B K (s)= 



1 

•so 

where r = 0.33, s = 37.46 GeV 2 , and 

r\ r 2 



Ff K (s) = 



1-4 (i--r) 2 ' 

m j v my 

where r x = 0.162, r 2 =0.173 and m x = 5.41 GeV. 



(27) 



(28) 
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Table 1: Leading order (LO) and next-to-leading order (NLO) coefficients af! , a P iv and a P iw entering into Eqs. d6lfTTTi. The NLO 
coefficients are the sum of the LO coefficients plus next-to-leading order vertex and penguin corrections. The superscript p is 
omitted for i = 1, 2, 3, 5, 7 and 9, the penguin corrections being zero for these cases. 
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f\ f \ A 1 'f\ 11/1 

-0.041 - zO.l 14 


a 3 
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-0.057 - /0.017 
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-0.062 - /0.004 
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0.0001 


0.0 + /0.0001 


< 


0.0007 


0.0008 + zO.O 


0.0007 


0.0007 + /0.0 








0.0007 


0.0008 + zO.O 
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-0.0094 


-0.0097 - /0.0003 


a" 
"lO 


-0.0009 


0.0005 + /0.0013 






-0.0009 


0.0006 + /0.001 




-0.0009 


0.0005 + /0.0013 






-0.0009 


0.0006 + /0.001 



Table 2: Parameters of production functions R"(E) and R^E) defined in Eq. (TBI for a- =3.506 GeV 
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a" 


r? (GeV- 1 ) 


col (G eV_2 ) 


a s 


•^(GeV- 1 ) 




1 


0.6731 


-0.2511 


1.5301 


0.3743 


-0.1090 


0.1008 


2 


0.6116 


0.0428 


1.5232 


0.7075 


-0.1029 


0.3256 


3 


1.2055 


0.3589 


3.1556 


1.0028 


+0.0979 


0.4653 



where cos0 g corresponds to the value of cos&h in Eq. ( 1321 ) 
with mi2 = W23. The helicity angle distribution dBr~ I d cos &h 
can be obtained from Eq. d29l by integration over the specific 
range of the effective mass m-a- 

Our aim is to describe the data of the Belle HI] and BaBar Q\ 
Collaborations in one common fit. The data chosen by us 
include the total branching fraction for the decay B ± — > 
[0(1020)^,0(1020) -h> K + K ], the averaged effective mass 
distributions dBr ± /dni23 for m.23 smaller than 1.8 GeV, and the 
averaged helicity angle distribution dB^ /d cos&h for 11123 < 
1.05 GeV. The distributions of the -> K+K K* events are 
obtained from the published data by subtraction of the back- 
ground components. The total number of data points for ten 
plots from both collaborations is equal to 175. The theoretical 
distributions are normalized to the total number of experimen- 
tal events corresponding to each data set. In our fit we used 
the averaged B ± -> [0(1020)/^, 0( 1020) -> K + K ] branch- 
ing fraction equal to (4.06 ± 0.34) • 10" 6 |7[]. There are four 
fitted parameters: ^, k, c and Np. The first three parameters 
are related to the S -wave decay amplitudes and the fourth one, 
Np, is the common f-wave normalization constant by which 
the amplitudes A p and A + p are multiplied. We have performed 
the fit to the 176 data points obtaining the total value of x 1 
equal to 343 and the following parameters: x - (6.44 + 0.44) 
GeV \ k = (3.51 ± 0.20) GeV, c = (0.084 ± 0.010) GeV 4 
and N P = 1.037 ± 0.014. For N P = 1 we obtain the aver- 
aged B ± [0(1020)^,0(1020) -> K + K ] branching fraction 
equal to 3.73T0 -6 which is within one standard deviation from 



3. Results 



Partial wave analysis of the decay amplitudes helps in the in- 
vestigation of the density distributions in the Dalitz diagrams. 
In Eqs. (12I3I > we have defined the S- and P- wave amplitudes to 
which the double differential B — > K^KiK^ branching frac- 
tion Br is related through the symmetrized amplitude A~ ym of 
Eq. (H2l): 



d 2 Br 



1 m 2 3\p\\\p2\ 

dm 2 3d cos&h E B 8(2tt) } mI 



U: vm (m 23 , ®h) 



(29) 



Here Tg is the total width of the B meson and the kaon mo- 
menta are: 

(30) 



lpil = 5 yf> 



4. 



\pt\ = VK - (mM + mK A i M l - (m23 - mK)2 \ (31) 

The helicity angle &h is kinematically related to the effective 
mass m\2 of the K7Kt system: 



1 



cos ti H = 



Wx\\p 2 \ 



-U M 2 B -m 2 2 ,+3m 2 K ) 



(32) 



Due to the symmetry of the Dalitz plot density under the ex- 
change of the kaons and K^, one can define the effective 
mass m23 distribution integrated over the ntu masses larger than 
m 23 : 

dBr f 1 d 2 Br 

-d cos&h, (33) 
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dm,2T,d cos &h 
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Figure 1: Moduli of kaon scalar non-strange and strange form 
factors (solid lines) obtained in our fit. The dashed and dotted 
lines represent the variation of their moduli when parameter k 
varies within its error band. 



the experimental value of (4.06 ± 0.34) • 10~ 6 . One sees that the 
absolute normalization of the P wave is very close to 1 which 
means that the decay amplitudes calculated in our model are 
adequate. 

Our value of x parametrizes a large range of K + K~ effec- 
tive mass up to 1.8 GeV and not just the region of /o(980). 
Therefore it cannot be directly compared with the value given in 
Ref. [4]. In addition, the estimate of^- given in Eq. (18) of Jj] 
involves the coupling constant of /o(980) to nn while here we 
have coupling to KK. Using gf KKlgf a nn - 4.2 from Ref. 
a very rough estimate similar to that given in Ref. ||4j]leads to 
X ~ 5.6 GeV 1 . The k parameter was not used in Ref. [4] where 
only the on-shell contributions to the form factors were taken 
into account. The value of k — 3.51 GeV 1 is reasonably larger 
than the typical KK mass considered. We have also done an 
analogous fit to the data using the three P-wave form factors 
based on the parameterization of Ref. (5|] obtaining similar val- 
ues of parameters as those written above, however with a higher 
X 2 value of 354. 

Fig. 1 shows the moduli of scalar form factors Y\ (523) and 
T s 2 (523) which determine the functional dependence of the S- 
wave amplitudes on the K + K~ effective mass. There are two 
prominent maxima of both form factors, one related to the 
/o(980) resonance and the second one forming cusps due to 
the opening of the third channel at 1400 MeV (in the present 
model responsible effectively for the production of four pions). 
Presence of /o(980) leads to the threshold enhancement of the 
S -wave amplitude. This effect can be directly studied in high 
statistics experiment with a very good effective K + K~ mass res- 
olution of about 1 MeV and should be seen only a few MeV 
above the threshold. 

In Fig. 2 the K + K~ effective mass distributions are shown for 
two mass ranges and for the data from the BaBar Collaboration. 
At low iuk+k- the spectrum is influenced by the P-wave ampli- 
tude and dominated by the <p( 1020) resonance. Above 1.05 GeV 
the S -wave amplitude is much more important than the P-wave 
one. According to our analysis which uses the approach of Refs 
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Figure 2: The K + K~ effective mass distributions from the fit 
to BaBar experimental data [2] in the 0(1020) range (a) and 
between 1.05 GeV and 1.8 GeV (b). Theoretical results are 
shown as solid line in (a) and as histogram in (b). 



I IJJ, |21j, the experimental maximum near 1.5 GeV can be at- 



tributed to the /o(1400 - 1460) found therein in solution A. We 
recall that in Ref. I2H1 the coupling constant of the /o(1400) 
decay to KK is much smaller than the corresponding coupling 
to nn. Let us notice that the model distribution depends on the 
sharp 4jt threshold located at 2 ■ m.3 = 1 .4 GeV which in real- 
ity should be smoothed out by the four-body pion interactions 
not taken into account in this quasi-two-body approximation. 
We have also studied the Belle Q_1J] K + K~ effective mass spec- 
tra and found that the quality of their description is similar to 
that shown in Fig. 2 for the BaBar data. Fig. 3 shows a more 
detailed comparison of the m^K- theoretical distributions with 
the Belle data [1], with events grouped in five ranges of 
which is the other combination of the K + K~ effective masses. 
One observes an overall general agreement of theoretical his- 
tograms with experiment, with some surplus of experimental 
events in Fig. 3e for the case of the highest slice of myi (larger 
than 20 GeV 2 ) where our model is not fully applicable due to 
the proximity of the Dalitz plot edge. 

Finally, in Fig. 4 we present the helicity angle distribution 
in the K + K mass range dominated by the 0(1020) resonance. 
Without the S -wave component of the decay amplitude the dis- 
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Figure 3: The K + K effective mass distributions from the fit 
to Belle experimental data JH (a) for m 2 n < 5 GeV 2 , (b) for 5 
GeV 2 < m\ 2 < 10 GeV 2 , (c) for 10 GeV 2 < m z n < 15 GeV 2 , 
(d) for 15 GeV 2 < m\ 2 < 20 GeV 2 and (e) for 20 GeV 2 < m\ v 
Theoretical results are shown as histograms. 



tribution should be symmetric with respect to cos ®h = 0. 
However, we observe an interference effect which distorts the 
distribution. This is a direct evidence of a non-zero part of the 
S -wave present even under the huge peak of the 0(1020) reso- 
nance. A theoretical integration of the S -wave contribution to 
the spectrum in the range from threshold till 1.05 GeV 

leads to about 12% relative branching fraction. It corresponds 
to the average branching fraction of 4.83T0 -7 which is in agree- 
ment with the experimental upper bound of 2.9 • 10~ 6 found in 
Ref. [1]. This agrees also with the BaBar estimate (9 + 6)% of 
the S -wave fraction in the region of masses between 1.013 and 
1.027 GeV Q. In the range of the K + K effective mass from 
1200 to 1800 MeV, which might be relevant for the X (1550) 
discussed in Ref. J^, the CP averaged branching fraction cor- 
responding to the S -wave is equal to 4.42 ■ 10~ 6 which is larger 
than the total contribution of the 0(1020) resonance. 

We have also studied the CP violation effects comparing the 
magnitudes of the decay amplitudes of the B and B + decays. 
While the moduli of the P-wave amplitudes for these charge 
conjugated decays are rather similar, the S -wave amplitudes 
behave differently indicating an important CP violation effect 
which depends on the mk+k- range. For the S -wave parameters 
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Figure 4: Helicity angle distribution for the Belle data [ 1] in the 
K + K effective mass up to 1 .05 GeV. The dashed line represents 
the S -wave contribution of our model, the dotted line - that of 
the P-wave, the dot-dashed - that of the interference term and 
the solid line corresponds to the sum of these contributions. 



written above, starting from the K + K~ threshold up to about 1 .4 
GeV, the modulus of the B + S -wave amplitude is larger than the 
corresponding modulus of the B amplitude. Then, above 1.4 
GeV, the B moduli become larger than the B + ones. Defining 
the CP asymmetry as 



lc/>(»*23) = ( 



dBr 



dBr + s fdBr 



dmn dmn ' v dmn dm-n 



dBr + s 



(34) 



one gets very large asymmetries if one takes into account solely 
the contribution of the S -wave. For example, A^ p (l GeV) = 
-0.51, A s cp ( 1.020 GeV) = -0.54, A S CP ( 1.25 GeV) = -0.95, 
and A^(1.50 GeV) = +0.59 (here the superscript S stands 
for the S- wave asymmetry). When the f-wave is included 
then the CP asymmetry is reduced to: Ac^(l GeV) = -0.25, 
A C p(1.020 GeV) = +0.029, A CP (1.25 GeV) = -0.85 and 
A C p(1.50 GeV ) = +0.495. Let us note a particularly small 
asymmetry in the range of the <p resonance, where the f-wave 
amplitude dominates, and an inversion of the Aqp sign above 
1 .4 GeV. Due to cancellations between the ranges of the nega- 
tive and positive asymmetries the resulting CP asymmetry av- 
eraged over the niK+K- range from threshold up to 1.8 GeV is 
rather small, equal to -0.05. The averaged branching fraction 
for the same mass range equals to 9.6 • 10~ 6 . It is worthwhile 
to add that the S -wave gives to it the dominant contribution of 
5.8 • 10~ 6 . 
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4. Conclusions 

We have studied final state interactions between kaons in the 
B ± — > ^K'K* decays. An overall general agreement with 
the Belle and BaBar data has been obtained. Our formalism 
is based on the QCD factorization supplemented with the in- 
clusion of the long distance K + K~ interactions. The latter are 
taken into account through the functional dependence of the 
scalar and vector form factors on the effective K + K~ masses. 
A unitary model is constructed for the scalar non-strange and 
strange form factors in which three scalar resonances /o(600), 
/o(980) and / (1400 - 1460) are naturally incorporated. The 
scalar resonance /o(980) leads to the threshold enhancement of 
the S -wave K + K~ amplitude. The K + K~ structure seen near 
1.5 GeV can be attributed to the third scalar resonance. A 
potentially large CP asymmetry is obtained in the mass spec- 
trum dominated by the S -wave. It originates from violent phase 
variations of the two kaon scalar form factors which affect the 
K + K~ effective mass dependence of the S -wave decay ampli- 
tudes. In general one can best study this effect away from the 
0(1020) peak. We have shown, however, that even under the <p 
maximum one observes nonnegligible helicity angle asymme- 
try. This effect originates from the interference between the S - 
and P- waves. 

Our approach presented here for the B ± — > K + K K ± de- 
cays can be extended to study the B° — > K + K K° S reactions 
for which results of the time-dependent Dalitz analyses have 
been recently published by the Babar [ 17] and Belle [18] Col- 
laborations . For further studies of the charged B decays new 
experimental data with better statistics are needed. Such data 
already exist! For example, the Belle Collaboration has now 
five times larger data sample than that used in their publication 
(jl analysed by us here. Future results from LHCb and from 
super-B factories would also be very useful. 
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